


X 



KEK-TH-964 
hep-ph/0406073 



Determination of 623 in Long-Baseline Neutrino Oscillation 
Experiments with Three-Flavor Mixing Effects 

Hisakazu Minakatcl and Masashi SonoyamctJ 

V:+". Department of Physics, Tokyo Metropolitan University, Hachioji, Tokyo 192-0397, Japan 

o: 

^ ■ Hiroaki Sugiyameo 

Q I Theory Group, KEK, Tsukuha, Ibaraki 305-0801, Japan 

Q ; (Dated: February 2, 2008) 

O ! Abstract 

We examine accuracy of ^23 determination in future long-baseline (LBL) v^ disappearance exper- 
ts ' iments in the three-flavor mixing scheme of neutrinos. Despite that the error of sin'^ 2^23 is indeed 
JC ■ a few% level at around the maximal mixing, we show that the error of physics variable s^g is large, 
C"^ ■ '^(•523)/'523 — 10-20%, depending upon regions of ^23- The errors are severely affected by the octant 
^ , degeneracy of ^23, and (^(slg) is largely amplified by the Jacobian factor relating these two variables 
f^ \ in a region near to the maximal mixing. The errors are also affected by the uncertainty due to 
^T ' unknown value of ^13; (^(slg) is doubled at off maximal in the second octant of ^23 where the effect is 
largest. To overcome this problem, we discuss combined analysis with t'g appearance measurement 



in LBL experiments, or with reactor measurement of ^13. For possible relevance of sub-leading 
effects even in the next-generation LBL experiments, we give a self-contained derivation of the 



Q \ survival probability to the next to leading order in sfg and Am|]^/A7TT-|^. 
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I. INTRODUCTION 

At several years since the pioneering discovery in atmospheric neutrino observation by 
Super-Kamiokande (SK) [l|, we now have compelhng evidences for neutrino oscillation. The 
direct confirmation of oscillation of atmospheric neutrinos comes from the K2K experiment 
using man-made neutrino beam [2] which confirmed its existence at 3.9cr CL J3|. The 
neutrino oscillation with the solar Am^ scale is beautifully observed by the KamLAND 
experiment j^ |5|. By observing deficit in reactor antineutrino fiux it pinned down the 
large-mixing-angle (LMA) solution of the solar neutrino problem based on the Mikheyev- 
Smirnov-Wolfenstein (MSW) mechanism [6|,|7|. It has thereby settled the long-standing solar 
neutrino problem J3] , whose origin was identified as due to the neutrino fiavor transformation 
by SNO |3] • Long awaited confirmation of the oscillatory behavior is now available by the 
three experiments, SK [lO|, K2K J3[, and KamLAND [5|. Thus, we already know roughly the 
structure of lepton fiavor mixing in the (1-2) and (2-3) sectors of the Maki-Nakagawa-Sakata 



(MNS) matrix [11 



We recognize that there are two small (or possibly vanishing) quantities in the lepton 
mixing matrix, ^13 and deviation of ^23 from 7r/4. The former is constrained mainly by the 
reactor experiments 12| to be sin^2^i3 < 0.15 at 90% CL by a global analysis [l3|, while 



the latter is bounded as < 1 — sin^ 26^23 < 0.08 at 90% CL by the recent reanalysis of 
atmospheric neutrino data by the SK group [ij]. Possible correlation between these two 
small quantities may imply hints for the yet unknown discipline that was used by nature to 
design the structure of lepton fiavor mixing. Some symmetries have been discussed by which 
the maximal ^23 is correlated with vanishing 6'i3 in the symmetry limit [15|. Furthermore, 
the question of how large is the deviation of 6*23 from 7r/4 is discussed from a variety of 
viewpoints |16]. Therefore, it is important to determine the value of not only 9i^ but also 
6'23 very precisely. 

How accurately ^23 can be measured experimentally? It appears that measurement of 
disappearance probability of z/^ at around the first oscillation maximum has highest sen- 
sitivity to sin^26'23.^ In the J-PARC -^ Super-Kamiokande (JPARC-SK) experiment, for 
example, sin^2^23 is expected to be determined to 1% accuracies [1^. Given such an enor- 
mous accuracy, it seems that there is not much to add because it looks like the best thinkable 
accuracy of mixing angle determination in neutrino experiments. Unfortunately, it is not 
quite true. We discuss three relevant and mutually related issues in this paper which prevent 
an accurate determination of 6'23, in particular, its small deviation from the maximal. The 
points we discuss are: 

• a large Jacobian associated with transformation of the variable from sin^ 2^23 to 5^3 
at around the maximal mixing 

• the octant ambiguity of ^23 (two solutions of ^23 larger or smaller than 7r/4) 

• effects of three-flavor mixing, in particular of non- vanishing 6*13 



^ The statement applies also to the other large mixing angle, 612, if v^, from reactors is used |l3|- In this 
case, the accuracy for sin 29i2 measurement translates into the one of sin 612 without suffering from the 
problems (apart from the third one) mentioned below. 



While the first point is "kineinatical" in nature, the second and the third points reflect 
the three-flavor structure of lepton mixing. Although each of these points was discussed 
in various occasions in the past, we believe that they are never discussed in a coherent 
fashion in the context of precision measurement of ^23 by long-baseline (LBL) experiments. 
The estimation of errors of ^23 determination is also done by taking account of three-flavor 
mixing effects [l9|,|20| with consistent results with ours. Yet, they offer neither understanding 



of the reasons behind the estimated large errors in the similar depth as we provide, nor ways 
of improvement by combining with e.g., reactor measurement. 

Importance of experimental determination of ^13 is now well recognized in the community 
by itself and as a door to exploration of leptonic CP violation. Strategies for determi ning 



the remaining oscillation parameters are developed with use of accelerator superbeam 21 1 
and reactor neutrinos [22j, which followed by a series of feasible experimental programs 18, 
23|, |2J, |25| • The accuracies of 613, determination in the LBL and the reactor experiments are 
expected to reach to sin^ 613 ^ 2.5 x 10^^. We feel that the strategy for similar accurate 
determination of 623 is relatively less well developed, and it is the purpose of this paper to 
trigger such efforts. 

By indicating that the three-flavor nature of the oscillation probability already manifests 
itself in the next generation LBL experiments, our discussion in this paper might have an 
implication to future neutrino oscillation research, suggesting a possible direction. The next 
to leading order corrections may be important due to the larger ratio of e = Am'^i/Am'^^ ~ 
1/30 due to a larger Am^i |5|], and possible large value of 6*13 close to the CHOOZ bound. 
Under such hope, we present in Appendix A a derivation of the full expression of survival 
probability -P(t'^ — > i^^) to the next to leading order in e and sf^, which is used in our analyses 
throughout this paper. See also J26| for derivation of such formulas under the similar spirit. 



II. WHY IS ACCURATE DETERMINATION OF ^23 SO DIFFICULT? 

In this section we present pedagogical description of the three reasons for difficulty of 
accurate determination of 6'23 that were mentioned in Sec. |l] The readers who are familiar 
to them are advised to skip this section and go directly to Sec. IIIII 



A. Problems of large Jacobian at around the maximal mixing 

The ultimate purpose of the whole experimental activities which aim at measuring neu- 
trino mixing parameters is to understand physics of lepton flavor mixing. For this purpose 
precise determination of the angles or the sine of angles are required. 

The z/^ disappearance measurement is sensitive to sin^ 26*23, not s^^, which is unfortunate, 
but is actually the case. It inherently makes accuracy of determination of 6*23 much worse 
than it is achieved for sin^ 2^23 at around the maximal mixing. Assuming that the true 
value of ^23 is vr/4, for example, an accurate measurement of sin^ 2^23 with error of 1% level 
translates into large error of ±10% level for S23. It is simply because the Jacobian for the 
transformation of the variables 



K43) _ dsls 



2 0/3 Ar.r.cOa^' ^ '' 



S{sm^ 2623) d sin^ 2^23 4 cos 2^23 
are large at around 6*23 = 7r/4, by which a tiny error obtained for sin^ 26*23, denoted as 



5(sin^26'23) in ((H), can be translated into a large error 6(823) ^^^ ■^23- I^ is what happens in 
Fig. El in Sec. Uni 

We emphasize that use of 5^3 , not sin^ 2^23, is inevitable in the error analysis based on 
three flavor mixing because the ly^ survival probability itself cannot be expressed in terms 
of sin^26'23 only once next to the leading order correction is taken into account. See (jHl) 
below and ()A15j) for the probabilities in vacuum and in matter, respectively. It is even more 
so if we combine Ue appearance measurement into our analysis, as will be done in Sec. IIV| 



because the dominant term in the Ue appearance probability depends upon 6*23 through s^ 



23- 



B. Octant ambiguity of ^23 

The dominant term in z/^ survival probability, and hence the disappearance measurement , 
is sensitive primarily to sin^ 2023- Therefore, there exist two solutions of ^23 for a given value 
of the survival probability -P(z^^ — > i^^), the octant ambiguity of 6^23.^ In Fig. [T] we present 
a contour of 90% CL allowed region in 523-5^3 space by 5 years measurement at JPARC-SK 
based on the analysis in Sec. II I II whose details will be explained in Appendix B. 

It should be noticed in Fig. ^ that there exist "clone" solutions, which arise as a con- 
sequence of the degeneracy of 6*23 solutions 29\. Namely, with an assumed true value 
S23 = 0.34, for example, the measurement allows two regions as solutions, one around the 
original value of S23 in the first octant and the other one around S23 = 0.66 which is the 
reflected point with respect to 6'23 = 7r/4 in the second octant. Of course, it does not produce 
additional experimental error of S23 as far as two regions are separated, though ambiguity 
exists of which region contains the true value of ^23- But once the two allowed regions start 
to merge (which takes place near the maximal 6*23), it does produce additional error of 5^3 
as exhibited in the spike in Fig. El The similar figure as Fi g. |T1 indicating ^23-^13 correlation 
and clone solution is used in a slightly different context in 



C. Three-flavor mixing effects in the disappearance measurement of 6 



23 



While the effect of three-flavor mixing gives relatively mild influence to the accuracy of 
6'23 measurement compared to the above two effects, it is nonetheless non-negligible. Let us 
explain the key point by taking a very simple case using the approximation of vacuum oscil- 
lation and small deviation from maximality of ^23- Though simplified, the model captures 
some relevant features possessed by the more realistic cases. 



^ It was recognized since sometime ago that there exists multiple solutions of mixing parameters, ^13, ^23 
and S, for a given set of measured oscillation probabilities. It is called the problem of parameter degeneracy 



and d tor a given set 01 
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I, and the octant ambiguity of 623 is only apart of it. The nature of the degeneracy 
may be characterized as the intrinsic degeneracy of 6*13 and 6 [23, which is duplicated by the unknown 
sign of Arn^i |23| and the octant ambiguity of 6*23 for a given sin 26'23 [2^ ■ The octant ambiguity may be 
the most difficult one to resolve am ong the eight-fold parameter degeneracies unless it is combined with 
either reactor 1221 or silver channel |33| measurement. 
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FIG. 1: 

ues given by nature are [(sls)^'^'^*' (•sfs)'^'^^*] = [0.34, 0.026] (red solid curve), [0.5, 0.026] (green 
dashed curve), and [0.66, 0.026] (blue dotted curve); Note that S23 = 0.34 and 0.66 correspond to 
sin^ 2^23 = 0.9 which is the lowest edge of the current allowed region and s^g = 0.026 corresponds 
to sin^ 2^13 = 0.1. The systematic errors are assumed to be (a^reh "^/^bg) = (10%, 20%). 



We use, throughout this paper, the standard notation of the MNS matrix j3, 
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where c,-,- and 
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difference of neutrinos is defined as Am 



The mass squared 

is the eigenvalue of the ith 

mass-eigenstate. Then, the disappearance probabihty 1 — -P(t'/i -^ i^fi) in vacuum is given 

by 
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and sin^jj, respectively. 
rn^ — rn^j where m. 
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In the right-hand side of Q, the first term represents major contribution which survives in 
the two- flavor or the single mass-scale dominance approximation J2a]- The second and the 
third terms are small corrections due to non- vanishing s\^ and e = Am\^/ Ain\^, respectively, 
which give effects of the three-flavor mixing. It is obvious that the second term vanishes for 
the maximal 6'23, while the third term vanishes at the oscillation maximum, Am\iL/2E = n. 
Suppose that we perform a hypothetical experiment at a monochromatic energy close 
(but not identical) to the oscillation maximum, sin{Am'^iL/2E) ^ 1. We assume that 



deviation from maximality of 623 is small, and parametrize it by t] defined as 77 = cos 2623 
(so that 1 — sin^ 2623 = r^^ -C 1). Using the approximation sin (Am'^^L/AE) = 1 we obtain 

P(^M ^ ^m) = Viv + Hs) + 2 4e sin [-^§-) ■ (4) 

Notice that the linear term in r^ in (Q distinguishes between 6*23 in the first and the second 
octants because f] changes sign depending upon which octant it lives. The equation (JH) 
clearly indicates that the larger value of 5^3 gives the smaller (larger) value of rj in first 
(second) octant of 6^23 for a given P(z/^ — > t'^). It is interesting to observe that the feature 
of 6'23-6'i3 correlation is in fact shared in Fig. ^ 

Therefore, what the disappearance measurement actually determines is not sin^ 2^23 but 
allowed regions in ^23-5^3 space. The feature leads to larger errors of ^23 in three-flavor 
analysis compared to those of two-flavor analysis. It is notable that the effect of 5^3 is 
slightly more significant in the second octant of 6'23. Note also the tendency of departure 
from the maximal ^23 due to ^13 is smaller (larger) at larger ^13 in the first (second) octant 
of 6'23. These properties are the key to fuller understanding of the features of the errors of 
6*23 to be discussed in the next section. T he p roblem of 6^13 dependence of accuracy of 6^23 
determination has been addressed, e.g., in J37|. 

Another notable feature in Q is that the solar mass scale correction can come into play 
a role. It can change the depth of the dip of P{i'fj, -^ i'^) and therefore may affect the value 
of 77. 



III. DISAPPEARANCE MEASUREMENT OF 6^23 AND THREE-FLAVOR MIX- 
ING EFFECT 

In this section, we fully analyze the sensitivity of ^23 by z/^ disappearance measurement 
by taking account of the three-flavor mixing effects. To make our discussion as concrete 
as possible, we take the particular experiment, the JPARC-SK experiment jlq|. While the 
Jacobian and the degeneracy effects are universal to any experiments, the three-favor effect 
is most prominent if the systematic uncertainties of the experiment are small enough to allow 
measurement of sin^ 2^23 to a few % level accuracies. We describe in detail the statistical 
method used in our analysis in Appendices B and C so that we can concentrate on physics 
discussion in the main text. Through the course of formulating our analysis procedure, we 
will try to reproduce the results of sensitivity analyses by the experimental group. (See 
Fig.El) 

We take 22.5 kt as the fiducial volume of SK and assume, throughout this paper, 5 
years running of neutrino mode v^ disappearance measurement, as expected in its phase I 
defined in LOI. We examine two cases of systematic errors, (cr^reh (^/^bg) = (10%, 20%) and 
(5%, 10%). The former numbers, which are based on their experience in K2K, can be too 
pessimistic but are quoted in LOI, while the latter may be a reasonable goal to be reached 
in phase I of the J- PARC neutrino project. To make the goal and to achieve even smaller 
errors the group plans to build an intermediate detector at ~ 2km from J-PARC We use, 
unless otherwise stated, Am|^ = 3 x 10~^eV^, Am^i = 7.3 x 10"^eV^ and s^2 = 0.315. We 
take the earth matter density p = 2.3 g ■ cm~^ based on the estimation quoted in [38|. 

We assume that Am\i > in our analysis. We have checked that the error (5(sin^ 2^23) 
slightly increases in the case of negative AttIj^ but only by about 0.001 in the allowed region 

6 



of 02^^^. In this paper, the quantity with superscript "best" always imphes the best fit value 
given by nature. We have also checked that our results do not change in any appreciable 
manner by taking Amli at the highest end of the LMA-I parameter, Amli = 10~^eV^. 
Therefore, we do not discuss further about the possible effects caused by fiipping the sign 
of Arrig]^ and by raising the solar Am\^. 

A. Accuracy of determination of S23 vs. that of sin^ 2^23 

We compare the accuracies of determination of sin^ 262Z and S23 to illuminate how the 
problem of Jacobian and the degeneracy affect the accuracy of 6'23 determination. We show 
in Fig. 121 the errors expected for determination of sin^ 2^23 as a function of sin^2^2r*. In 
Fig. I^t' and Fig. |2b, the systematic errors are taken as (cr^rei, cT/^bg) = (10%, 20%) and (5%, 
10%), respectively. The figure includes not only the result of the three-fiavor analysis with 
the particular nonzero input value of ^13, sin^ 2^^!*^* = 0.1, but also that of the two-fiavor 
analysis assuming ^13 = 9\l^^ = 0. We note that the results depend on input values of 6^^^^ 
but only mildly; If we take sin^ 29\^^^ = in the three fiavor analysis, the uncertainty of 
sin^ 26^23 decreases only by ~5%. 

We call special attention to the distinction between the lower solid lines, and the upper 
dashed and dotted lines in the three-fiavor results in Fig. |21 Namely, the red thin-solid and 
the pink thick-solid line are obtained under the assumption that 6'23 is in the first octant, 
< ^23 ^ 7r/4, by which the octant degeneracy is switched off by hand. This procedure is 
adopted, for example, in LOI of the JPARC-SK experiment in their sensitivity estimate |l8[. 
Whereas for the green dashed and the blue-dotted lines, the full region < 6*23 < 7r/2 is 
allowed as output, and the input value of 6*23 is taken in the first (second) octant of 6*23 for the 
green dashed (blue dotted) lines. We recognize that our ignorance of which octant 6^23 lives 
makes the error of sin^ 26'23 almost twice (~ 80% increase) in the optimistic (pessimistic) 
case of systematic errors. 

In Fig. 121 the errors expected for determination of S23 are plotted, with systematic errors 
(o'^rei, ct^bg) = (10%, 20%) and (5%, 10%) in Fig. 12^- and Fig. |2b, respectively. The readers 
may be surprised by the difference between Fig. |21 and Fig. |21 Although they are based on 
exactly the same data set, the errors S{s23)/s23 in Fig. El are as large as 14%-20% in the 
"near to the maximal" region and the shape is quite odd, which should be compared to 
the smooth behavior of the curves in Fig. |21 Hereafter, we use the notations "near to the 
maximal" and "off the maximal" regions which mean the regions of ^23, 1-523 "~ 0.5 1 — 0-08 
and 0.08 < |s23 ~ 0.5| < 0.16, respectively. There are two main reasons for such a marked 
difference. 

(a) Jacobian effect 

As explained in Sec. Ill Al the Jacobian (Q) makes the uncertainty of S23 much larger than 
that of sin^ 26'23 near the maximal mixing. It explains the feature that the uncertainty is 
large in region 0.45 < S23 < 0.55, but it does not explain the double-peaked structure in 
Fig. El 

(b) Degeneracy effect 

The problem of degeneracy related with ^23 requires a special comment. As one can 
see in Fig. HI there is a clone solution for a given input value of 6^23. We do not include 
the separated clone region in defining the uncertainties in SI3. We think it appropriate not 
to include the separate region because it is due to the problem of parameter degeneracy 
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FIG. 2: The expected accuracies of determination of sin^ 2^23 at 90% CL obtained by the three- 
flavor analysis are compared with that of our two-flavor analysis (^13 = O^^,^^ = 0) which is indicated 
by pink thick-solid line. The red thin-solid line is the result obtained in the analysis assuming that 
the ^23 is in the first octant, < ^23 < 7r/4. The green dashed and blue dotted lines are for input 
values of ^23^* i^ the first and the second octants, respectively, while allowing the output ^23 in tbe 
full range, < ^23 < '?'"/2. Fig. |2K and [213 are with systematic errors (a^reij ct^bg) = (10%, 20%) 
and (5%, 10%), respectively. 



and is a separate issue from the uncertainty attached to experimental S23 determination. 



Nevertheless, it affects the results in Fig. El because the genuine and the clone regions merge 
near the maximal mixing. Once the true and the fake regions merge there is no way to 
separate them. It creates a jump of S^s^^) at s^g ~ 0.45 and ~ 0.55 when the true value of 
S23 approaches to 0.5 from above and from below, respectively, leading to the double-peaked 
structure in Fig. El^ 



^ Strictly speaking, it may not to be appropriate to connect the points of S{s2^) before and after the merge 



0.14 

0.12 

0.1 

-^ 0.08 

CMCVJ 

'<' 0.06 

0.04 

0.02 





(MC\J 



0.3 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




0.3 



Am''=3.0x10"''(eV^),o, 



jirel=^0%, cj^Bi3=20% 



(a) 



3-flavor, sin^2e^f =0.1 
3-flavor, sin^2e^f' =0 



2-flavor 




0.4 



0.5 

/_2 ^best 



0.6 



0.7 



Am^=3.0x10"^(eV^), aj,rer5%, aj,BQ=10% 



(b) 



3-flavor, sin^2e^f =0.1 

3-flavor, sin^2e^3® =0 

2-flavor 




0.4 



0.5 

/_2 ^best 



0.6 



0.7 



FIG. 3: 



The expected accuracies of determination of s^ 



23 



at 90% CL obtained by the three-flavor 



analysis are compared with that of the two-flavor analysis (^13 = Oi^, = 0) which is indicated 
by pink thick-solid line. The red thin-solid and green dashed lines are for input values of ^13 as 
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f^MBc) = (10%, 20%) and 



= 0.0, respectively. Fig. |2K andlHt" are with systematic errors (cr^rei, 
10%), respectively. 



Notice that the sensitivity curves in Fig. |21 obtained not only with but also without the 
restriction < ^23 < 7r/4 do not have the similar jump. It is because in the analysis of error 
defined in terms of the variable sin^ 2^23 the true and the fake regions start to overlap at 
even larger deviation from the maximal value of ^23, leading to a smooth behavior. It stems 
from the fact that the variable sin^ 2^23 cannot distinguish between the first and the second 
octants, through which the confusion due to ^23 ambiguity is maximally enhanced. It also 
implies that use of the variable sin^ 2^23 in defining errors leads to larger errors in regions 
where 5 (533) jumps down in Fig. El We emphasize again that one has to use the variable 



in Fig. 131 because the rise of (^(523) is a discrete jump 



S23 to represent uncertainty in ^23 determination because, in addition to the one mentioned 
just above, the angle is defined in the entire quadrant and the obtained uncertainties are 
different in the first and in the second octants. 

From Fig. |21we also observe the following features. 

(i) The infiuence of the systematic error on the uncertainty in S23 determination is rather 
mild. With a factor of 2 smaller systematic errors, S{s2s) decreases only by ~15% in the 
near to the maximal region and by ^20% in off the maximal region. 

(ii) The three-flavor (nonzero 613) effect exists in the entire range of S23 within the SK allowed 
region. Its magnitude differs depending upon which region of s^^ and which octants of 623 
is talked about. Roughly speaking, the three-flavor effect makes (^(523) increase by ~ 0.015 
(20%-25%) in the near to the maximal region and ~ 0.008 (^ 0.017) in off the maximal 
region in the flrst (second) octant of 6*23. Note that the uncertainty is almost doubled in the 
second octant. 

(iii) It is indicated in Fig. 01 that the shape of uncertainty curve is distorted in the case of 
input value 8110^29^^^^ = 0.1; The "Mexican hat" moves slightly to rightward, and curiously, 
better sensitivities are obtained in some regions. Peak in the uncertainty curves signals 
that merging of the degenerate solutions occurs, and its exact location is determined by a 
delicate balance between the parameters. Therefore, the values of uncertainties around the 
peak region cannot be trusted at a face value. However, we have checked that the uncertainty 
at {823)^'^^^ = 0.5 does not suffer from the degeneracy problem in a severe way, as is shown 
in Fig. 121 Hence, we have restricted ourselves into a roughly estimated uncertainties in the 
above discussions. The uncertainties are not so sensitive to the input value of ^13 apart from 
this sensitivity to merging degenerate solutions. Especially, the uncertainty for (^23)'^^'^' = 0.5 
is almost independent of (^^3)^'''^* as we can understand with (JH). 

We have also examined the case of ^m'^i = 2.5 x lO^^eV^ to make comparison with 
that of A?7i|^ = 3 X 10~^eV^. It may worth to examine because it may mimic the situation 
that the JPARC-SK experiment run at slightly higher energies than the vacuum oscillation 
maximum. While the qualitative features noted above also hold in this case there are some 
notable changes: 

(iv) With Am|^ = 2.5 x 10~^eV^, the uncertainty in S23 determination becomes larger by 
~15% in the near to the maximal region and by ~20% in the off the maximal region. The 
feature is essentially the same for both the optimistic and the pessimistic systematic errors, 
(a^rei, (t^bg) = (5%, 10%) and (10%, 20%). 

IV. TO WHAT EXTENT DOES LBL Ue APPEARANCE OR REACTOR MEA- 
SUREMENT HELP? 

A better sensitivity for 6*23 can be achieved by combining some other experiments that 
are sensitive to ^13 with disappearance measurement, though it cannot solve the problems of 
the large Jacobian and the degeneracy; An improvement can be expected to the extent that 
the loss of sensitivity of 623 is caused by the uncertainty of 6^13. There are two possibilities 
for the supplementary experiments: z/e appearance measurement in LBL experiments and 
^13 measurement in reactor experiments. Fortunately, the z/^ disappearance measurement 
is always accompanied by the simultaneous z/g appearance measurement of Pii^i^ — »■ z/g) in 
most of the LBL experiments. Therefore, it is natural to think about combining them flrst 
to determine 6'23 more precisely. But if a reactor measurement of 6*13 is carried out it, it may 
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FIG. 4: The expected accuracies of ,523 determination at 90% CL when combined with i/g ap- 
pearance (green dashed hne) or reactor measurement of ^13 (blue dotted hne) . They indicate the 
improvement over the one reachable by three-flavor analysis of v^ disappearance measurement only, 
assuming sin^ 2^13 < 0.15, which is indicated by red thin-solid line. The true value of sin^ 2^13 is 
assumed to be 0.1 and 0.02 in Fig. ^ andlljj, respectively. The result of two-flavor analysis (pink 
thick-solid line), which is given for comparison, is almost degenerate with green dashed and blue 
dotted lines. 

give a better determination of ^13 than LBL i/^ appearance experiments because the latter 
suffers from the additional uncertainty due to 5 31 . 

Here, we present only the result for 6^^^^ = 0. We have checked that when the input 
value 6^'^^^ is varied in the range — 7r/2 < 6^^^^ < tt/2 the variation of ^(323) is within ±5% 
relatively to the value for 6^^^^ = 0. Note that the Ax^ must be minimized with respect to 
5 because we can not know the value of 6 from Ue appearance only. Again, we postpone the 
detailed description of our statistical procedure in Appendix C. 

The obtained results of how the sensitivity improves with LBL Ug appearance and reactor 
measurement are presented in Fig. I^for two values of 0\^^^, sin^ 2^^!*^* = 0.1 (Fig. 0^) and 
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0.02 (Fig. EJd). We examine the cases with two true values of sin^ 2^13 because the accuracy 
of sin^ 2^13 measurement in appearance experiments depends on its true value. On the other 
hand, the error of sin^ 2^13 in the reactor-neutrino disappearance measurement is nearly 
independent of the true value of ^13, 6 (sin^ 26'i3) ~ 0.02 at 90% CL, as was shown in our 



previous analyses [22,|47|. We utilize this feature to simplify our formulation of the statistical 
procedure with reactor measurement in Appendix C 

It is shown in Fig.|3]that the accuracy of 623 measurement is indeed improved thanks to the 
help by z/e appearance experiment or reactor measurement of 613. The reason why the better 
sensitivity of 623 is obtained at smaller sin^ 26*5^3^* in the case with appearance measurement is 
that the error of sin^ 2^13 determination increases for greater values of sin^ 29\f^^ due to larger 
CTereii^e^^^)'^ in ()C2|) . In fact, the case presented in Fig. |1|d with sin^ 26\^^^ = 0.02 represents 
the best sensitivity attainable by combined appearance-disappearance measurement. Notice 
that the error of 623 cannot be smaller than that of two-flavor fitting (^13 = 0) because it 
is the limit of no uncertainty in 6*13. Moreover, appearance measurement itself is not very 
sensitive to 623 and hence it plays only a supplementary role by constraining 613. 

It is also shown in Fig. |3]that with reactor measurement of 613 the error of 623 determi- 
nation is reduced to the one expected in the two-fiavor fitting with 6*13 = independent of 
the values of 6\^^^. Thus, the reactor measurement of 613 can do a better job in improv- 
ing accuracy of 623 determination in relatively large O^^'^^. It is due to the feature that the 
uncertainty in ^13 determination in reactor experiments is insensitive to the true values of 
^13, and that it is a pure measurement of ^13. In the combined analysis with appearance 
measurement, by contrast, min^Axg (see Appendix C for definition) vanishes not only at 
the best fit point but also in region around it due to the minimization in terms of 6, and no 
improvement of Axt is obtained in the region. 

We did not attempt to further combine Ue appearance measurement, because the sensi- 
tivity in two-fiavor analysis, which can be regarded as the limiting sensitivity attainable by 
the present method, is already achieved with reactor measurement of 613. 



V. CONCLUDING REMARKS 

In the next generation LBL neutrino oscillation experiments, a precise determination of 
sin^ 26^23 to a few % level is expected. In this paper, we have pointed out that, despite such 
marvelous sensitivities for sin^ 26'23, the reachable accuracy of determination of 5^3 is only 
at a level of S{s23)/sl^ = 10%-20% depending upon the true value of 6*23. 

We have identified the principal reasons for the disparity between the errors of sin^ 26'23 
and S23 as due to the Jacobian and the degeneracy effects. We also noticed that measurement 
of 6'23, when it goes down to such high accuracy, is started to be affected by the higher order 
corrections in sf^, in the oscillation probability. We have shown that by doing detailed 
quantitative analyses such three-flavor effects produce, typically, an increase of uncertainty 
in S23 determination by 0.008-0.017; 6(323) i^ larger in the second octant of 6'23 than that in 
the flrst octant. 

Then, we have shown that by combining either reactor measurement of 613, or Ue appear- 
ance measurement in LBL experiment the accuracies of 623 determination can be improved 
to the level expected by simple two-flavor fltting {613 = 0) of the measurement, the limiting 
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sensitivity attainable by the disappearance measurement.^ 

We have to note that our simulation of data and the estimation of uncertainties were 
done with OA 2 degree beam, which may be slightly higher in peak energy than the one 
appropriate for Am'^i = 2-2.5 x 10~^eV^. It is worthwhile to repeat the analysis with more 
realistic OA 2.5 degree beam. 

Finally, we want to make some remarks on importance of precise determination of the 
lepton flavor mixing parameters. We want to stress two viewpoints, one related to physics 
and the other to future neutrino experiments. 

(1) Why is precise measurement of the lepton flavor mixing parameters so important? 
One of the significant features is that there are two large and small mixing angles. It is 

suggested that near maximal 623 may imply a symmetry and 623 is maximal and 613 vanishes 
in the symmetry limit !l5|. If it is indeed the case detecting small deviation of 623 from the 
maximal is as important as measuring 613. Confirmation of this possibility would require 
a few % level determination of S23, given already severe constraint by the global analysis, 
■5i3 < 0.047 at 3cr CL [l3|. If it turned out that ^13 is much smaller, it would necessitate 
even more accurate determination of 623. If some of the lepton and quark mixing angles are 
indeed "complementary" with each other |32l], there will be immense requests for precise 
determination of lepton mixing angles because the quark mixing angles are now measured 
quite accurately. This possibility may require even more precise determination of 623 beyond 
a level expectable by current technology. 

Keeping in mind the possibility that we finally uncover the real theory of flavor mixing 
in the future it is desirable, whenever opportunity exists, to pursue determination of mixing 
parameters as precisely as possible. Notice that there is a way to measure 612 to an accuracy 
of ~ 2% comparable to that of the Cabibbo angle jl7|. 

(2) How important or crucial is the precise determination of 623 for future neutrino oscillation 
experiments? 

Most probably the field where precise determination of lepton mixing parameters is really 
necessary is the detection of leptonic CP violation, one of the ultimate goal of future neutrino 
experiments. It is because CP violation is a tiny effect, being suppressed by the two small 
numbers, one the Jarlskog factor ~ 0.04 sin 5 and the other e = Am2i/ ^m'^i ~ 0.033. Hence, 
it is crucial to minimize uncertainties of other mixing parameters (as well as neutrino cross 
sections on the relevant nuclei) to achieve unambiguous detection of CP violation. 

In this regard, diminishing the uncertainty in 623 measurement is probably one of the 
most important. To obtain the feeling let us suppose for simplicity that the experiment 
has relatively short baseline so that the vacuum oscillation approximation applies, which 
may be the case for the JPARC-Hyper-Kamiokande project 40]. Let us assume that the 
CP asymmetry defined hj A = {P^i^) — P{i'))/{P{v) + P{v)) gives a good measure for CP 
violation, where -P(j^) and P^i') denote appearance probabilities P{i'^ -^ I'e) and P{Ufj_ -^ z/e), 
respectively. To simplify the equations we tune the energy at the first oscillation maximum, 
/\m\'^L/2E = vr, as planned in the J- PARC experiment. Under these approximations, the 



^ In fact, accelerator ^13 measurement is affected by the disappearance experiment through S23 and/or 
Ato|]^. Therefore, we may need simultaneous determination of ^23 and ^13 for the most stringent bounds. 
While we did not go into deep into the detail on how such analysis for simultaneous determination of two 
angles can be done, such a proper formulation is necessary and should be ready when the JPARC-SK 
experiment starts. 
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CP asymmetry is given by 



P{v)-P{v) sm2e^2Cote23 fAml,L\ . 
A = ^, , — — — = — smd. (5) 

At the moment, of course, the uncertainty of S13 is a dominating error because it has not 
been measured yet. But it is expected that it can be measured either by LBL or reactor 
experiments with errors ^ 10-25% if sin^ 2^13 > 0.04. The error produced by cot ^23 factor 
in (0) is non-neghgible compared to this error. The uncertainties estimated in Sec. IIIII can 
be translated into ~ 5-15% error in cot6'23, and hence it is large. Here, it is worth to note 
that where the errors are relatively small the splitting of cot 623 between two degenerate 
solutions of 623 is large. Therefore, the uncertainty in the CP asymmetry is even larger, 
20-30%, in off the maximal region if the octant degeneracy of 623 is not resolved. Most 
likely, the error of S13 is significantly reduced at the time of CP measurement because of 
the possibility that it is precisely determined by in situ LBL measurement employing both 
neutrino and antineutrino beams at the oscillation maximum (or more precisely the shrunk 
ellipse limit |3l|). Therefore, the uncertainty of 623 including the degeneracy effect may be 
the dominating error for CP measurement. 

Precise determination of Am|^, though not explicitly addressed in this paper, is an an- 
other important issue. The error is expected to be ~ 10% in the JPARC-SK experiment. Its 
precise determination is also essential in a variety of contexts ranging search for CP violation 
and improving bounds on mixing parameters by double beta decay experiments. 

In conclusion, we stress the importance of finding alternative methods which are sensitive 
to S23 not to sin 2^23 for accurate measurement of 6^3. The obvious candidates include the 
reactor-LBL combined measurement discussed in [2^, and high-statistics observation of 
atmospheric neutrinos |4lL 1421 . 
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APPENDIX A: DERIVATION OF DISAPPEARANCE PROBABILITY P(u„ 



TO THE NEXT TO LEADING ORDER IN Aml^/Aml^ AND s 



In this Appendix, we give a self-contained discussion for deriving the expression of dis- 
appearance probabihty P{i^fj, -^ y^) which is vahd to the next to leading order in e and s^g. 
We use the method developed by Kimura, Takamura, and Yokomakura (KTY) 43] . 

The evolution equation of neutrinos can be written in the flavor eigenstate as 



d 
ax 



2E 



HaP^jS 



{a,p = e,fj.,T), 



(Al) 



where the Hamiltonian is given by 



H = U 
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(A2) 



whose first term will be denoted as H'""''^ hereafter and Ar/i?, 



jt 



m^ 



mf. (Hence, Amli = 
by definition.) In ()A2|) . a = 2\/2G pN f,{x)E denotes the index of refraction of neutrinos 
in medium of electron number density Ne{x), where Gp is the Fermi constant and E is 
the neutrino energy J6|. Despite that Ne{x) may depend upon locations along the neutrino 
trajectory, we use constant density approximation throughout this paper. The MNS matrix 
U relates the flavor and the vacuum mass eigenstates as 



l^n 



{U)ail'i 



where i runs over 1-3. 

We now define the mass eigenstate in matter i/™ by using transformation 



(A3) 



(A4) 



where V is the unitary matrix which diagonalize the Hamiltonian with scaled eigenvalues A 
as V'^HV = Hdiag = diag{aXi, aX2, a A3). We first obtain the expressions of the eigenvalues 
of the Hamiltonian ()A2|) . They are determined by the equation det[H — Xal] = which is 
the cubic equation for the scaled eigenvalue A: 

A^ - {621 + 631 + 1)X^ + [(c?2C?3 + sl^)62i + 4^31 + 621S31] X - 01^01^621631 = 0. (A5) 

Notice that everything is scaled by a and 621 and ^31 denote the scaled squared mass differ- 
ences, 



_ Am|i 
021 = ■ 
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31 



Amii 



(A6) 



To first order in ^21 and sf^ the solutions of the equation are given under the convention 
that Al < A2 < A3 by 
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Notice that we are in the intermediate energy region between high-energy (atmospheric) 
and low-energy (solar) level crossings. If we sit on the high-energy region above two level 
crossings, the expressions for A2 and A3 must be interchanged. We have checked that the same 
analytic formulas are obtained, when expressed in terms of observable physical quantities, 
even if we work in the high-energy region. 

We follow the KTY method for deriving P(z/^ -^ z/^) and write down the equations 



H, 



^ee^TT 






TTvac 



Tvac 

TT 



Tjvac Tjvac 



[HIT + o)Hl 
They give relationships between mixing matrix in vacuum and in matter as 



E^^i^-i' = E 



I ^ /j« I 
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(A8) 
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Solving ()A9|) for |V^2p under the constraint of unitarity 

E 1^-1^ = 1 



(A9) 



(AlO) 



we obtain 



W -l^ 



A,iA 



ki 



[XiP + q- Xi{Xj + Afc)] 



where A^j = Xj — Aj and {i,j, k) are cychc. 

The disappearance probability P{i'f_i — > t'^) is given by 



P(iy„ 



where Lj, is given by 
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Xi{Xj + Afc) + Aj(Aj + Afc) + AiAj(Aj + Xk){Xj + A^). 
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Notice that p and q are expressed by the vacuum mixing parameters through ()A9|) . If 
we insert the exact form of the eigenvalues Aj in ()A13|1 we obtain the exact expression of 
the disappearance probability Pii'fj, —>■ i^fj,), which should be identical with 1 — P{i'^ —>■ 
Ve) - P{y^ -^ Ur) computed by KTY 0. 

In this paper, we restrict ourselves into the perturbative expansion and keep terms to 
order 621, sfg, and S21S13 for CP 6 sensitive terms. Using the perturbative result of the 



eigenvalues ()A7|) we obtain, for example, 
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Collecting these formulas, we obtain the expression of the disappearance probability to the 
next to leading order as^ 
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where J^ = C12S12C23S23S13. The leading order term appears only in D^i (the first term). 
Notice that once we employ perturbation expansion by ^21 = AttiIi/o we cannot take the 
vacuum oscillation limit a — > in ()A16|) and ()A17|) . 

We also note a pathology in doing perturbative computation of the oscillation probabili- 
ties. In energy region around dip at the oscillation maxima, the leading term in the survival 
probability vanishes and the perturbative expansion is not well defined. In fact, the pertur- 
bative formula for P(z/^ — > I'fj,) in ()A15|) becomes slightly negative in the energy region. To 
obtain better approximation we may have to sum up small terms that are negligible except 
for this region. In the context of the analysis in this paper, however, we have explicitly 
checked that the pathological feature does not affect our estimation of the sensitivity. 



APPENDIX B: METHOD FOR STATISTICAL ANALYSIS OF DISAPPEAR- 
ANCE EXPERIMENT 

We define our statistical method to quantify our sensitivity analysis for determination of 



6*23. We use the following standard form of Axi'- 



AxJ(s23,Am2^,si3) = Y^ 



{N,^ 



Arbcst 



r 



N^"^ + B,^ + {a,,eiKn' + {cj,^oB,^Y 



(Bl) 



where N^i denotes the number of v^ charged- current quasi-elastic events of the reaction 
v^ + N —^ ^~ + N' m iih. energy bin which is to be tested against the artificially created 



^ While we were writing the manuscript of this paper, we have informed |44J that a similar perturbative 
formula was obtained by Akhmedov et al. 26]. We thank Tommy Ohlsson for kindly informing us the 
paper prior submission to the Archive and for verifying the consistency between their and our results. 
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"experimental event number" whose best fit is N^^f^^. B^i is the number of background events 
calculated in the similar way. We use 4 energy bins of width 0.2GeV in the energy range 
E„ = 0.4-1.2GeV. a^rei and ct^bg represent the systematic errors associated with signal and 
background events, respectively.^ 

To calculate the numbers of signal events we convolute the survival probability -P(t'^ — >■ 
z/^) with the neutrino flux times cross sections given in J45|. We use for Pii^^ — > u^) 
the expression valid to next to leading order in e and s^g, ()A15|) with ()A16|) and ()A17|) . 
We have explicitly checked that (5-dependent terms gives a minor effect in our analysis 
because they are of the order of esis which are higher order than terms we are keeping 
(assuming e and sfg are comparable in size). We assume 5 years running of neutrino mode 
Ufj, disappearance measurement. The estimation of background is done by using energy 
distribution of background events after cut kindly provided by the JPARC-SK group ^] . 
We have cross-checked our procedure against the numbers of events after cut. 

The procedure of calculating the error of 623 determination is as follows. We take an input 
value of 6^13 as experimental best fit value (or nature's choice), and obtain 90% CL allowed 
region in s'^^-Am^-^^, or sin^ 26'23-Am3]^ plane. To obtain the allowed region, we project three- 
dimensional manifold of Axt = constant onto the above two-dimensional plane by imposing 
the CHOOZ constraint 5^3 < 0.04 (sin^ 2^13 < 0.15) at 90% CL. We rely on the analysis of 
2 degrees of freedom because of the resulting two-dimensional plane. Notice that while we 
take a particular value of ^13 as an input a wide unlimited region is allowed for output value 
of ^13, as shown in Fig.^ The feature arises because the disappearance measurement is poor 
at restricting 6*13. Then, (5(sin^ 26^23) and S{sl^) plotted respectively in Fig. |21and Fig. |21are 
half width of the projected allowed region onto the sin^ 2^23 and S23 axis for a value of Am"^^. 
We examine cases of two input values of 6^13, sin^ 29^1^''^ = 0.0 and 0.1 ((^^3)^*'''* = 0.026). 

In the JPARC-SK's LOI jl8| they describe the procedure of determination of 623 by the z/^ 
disappearance measurement at the neutrino energy that corresponds to the first oscillation 
maximum. In Fig. 11 of [l8| (original version) they present, based on analysis assuming 
the two-flavor mixing, the accuracy of determination of sin^ 26'23 and Ani^i as a function of 
Arrtg^. To check the consistency of our procedure of computing event numbers as well as 
the method for statistical analysis we try to reproduce the former figure and the result of 
this analysis is given in Fig. El Fig. Et and Fig. Eb are for the cases sin^ 26'2|'^* = 1.0 and 0.9, 
respectively, and we have used the two-flavor analysis. We restrict ourselves to the plot of 
uncertainty of sin^ 26*23 because it is of importance for our later analyses. We use cr^rei = 0% 
and cr^BG = 20%, and take the earth matter density as p = 2.8 g • cm~^ as in LOI.^ Our 
results reproduce well the error of sin^ 2^23 presented in LOI. The agreement is particularly 
good for sin^ 2^2^* = 1.0 and is reasonably good for sin^ 2^23^* ~ 0-9 i^ ^^^ interesting region 
of Aml^. 



^ cr^rci denotes the relative normalization error between the numbers at front and far detectors, which 
corresponds to V2 times the uncorrelated error (e.g., of flux) between those detectors. We assume, for 
simplicity, that the baseline length of the front detector is short enough to make the numbers of events 
independent of oscillation parameters and neglect the background as well as the statistical error for the 
numbers of events in the detector. 

^ Despite the impression one has from the description in the LOI, they assume a^rei = in their plot in 
Fig. 11. We have explicitly confirmed this point by communications with the experimental group. We 
remark that it might not affect the results too much due to small number of events at the dip, S{N) = aN. 
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ArPgi^ vs 5(sin^2e23), sin^2e^§^'=1 .0, 0^^61=0%' 0^60=20% 
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FIG. 5: Comparison of our estimate of sensitivity for sin^ 2^23 (red solid line) with that of the 
JPARC-SK group given in LOI (blue dotted line). Both estimates use the two-flavor mixing and 
the same systematic errors. The input value of 023 is taken as maximal, sin^ 2023*^* = 1, in Fig. EK 
and as sin^ 26l^f * = 0.9 in Fig. Eh- 

APPENDIX C: COMBINED ANALYSIS WITH THE LBL APPEARANCE OR 
THE REACTOR EXPERIMENTS 



We define Ax^ for the combined analysis with the LBL Ug appearance experiment as 



^X^e(-S23,Am3i,si3) = Ax^(s23, Amg^, S13) + mm Axe(si3, 5, S23) 



(CI) 



^xl{Si3,6,S23) = 



{Ne-N^""'')' 



ATbcst + B, + aLi(iVe'-*)2 + 0-23^52 • 



(C2) 



The numbers of signal and background events are denoted in ()C2|) as A^e and B^,, respectively, 
which are defined without binning for Ue appearance. In the analyses in this section, the 
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systematic errors in disappearance and appearance measurement are assumed to be the 
same, cTerei = o'fj.rei = 5% and (TeBG = o"mBG = 10%. In the definition of Axg, we fix for 
simphcity 623 and Amg^ at their input values ^23*^* and {Arn^i)^^, respectively. For more 
details on our treatment of the appearance experiment, see J47|. The allowed region is 
obtained by projecting three-dimensional manifold of A^^^ = constant onto the s'^^-'^tn'^^ 
plane, the similar procedure as we defined in Appendix B for the analysis in Sec. 11111 

Here are some technical comments on how we simulate the number of events after exper- 
imental cut. Reconstruction of neutrino energy can be done for quasi-elastic events but not 
for inelastic events unless reaction products deposit all the energies and they are detected in 
the water Cherenkov detector. In our analysis, the relation between reconstructed neutrino 
energy Ej-^c of background events and true energy E„ is assumed to be 

Erec = E^~ O.SGeV (C3) 

by using the comparison between Monte Carlo simulation done by the JPARC-SK group. 
See Fig. 2 in 0. 

For signal events, Er^r = E^ is assumed. Then, from a spectrum Nac calculated by 
JPARC-SK group |46[ with a certain set Xq of parameter values, the spectrum after cut for 
any values of parameters x is estimated on a bin-by-bin basis as 

Nk ix) 
Nacix) = -^ N^xo), (C4) 

where Nbc represents the calculated number of events before cut. In our analysis, this 
prescription is used for all number of events in JPARC-SK. 

We define Ax^ for the combined analysis with reactor ^13 experiment via the following 
considerations: The atmospheric Am^^ will be determined by the location of the dip in 
the energy distribution of muons in disappearance measurement corresponding to the first 
oscillation maximum. Then, the reactor experiments can determine ^13 independent of the 
other oscillation parameters as well as the matter effect as discussed in |2^|. Therefore, it 
acts like giving additional constraint purely on ^13. The appropriate Ax^ then takes a very 
simple one, 

Axlr{s23, Am^g, S13) = Axl{s23, Aml^, 3^3) + Ax'(si3). (C5) 

For Axl of reactor experiment we use an effective parametrization 



^r 



rather than dealing with a particular experimental setup explicitly, where sin^ 2^5'3'^* is the 
true value of sin^ 2^13 assumed in the analysis. Note that the uncertainty in sin^ 26'i3 in the 
reactor-neutrino disappearance measurement is almost independent of true value of ^13 as 
was shown in our previous analyses [22,143]. This feature is in contrast with the one in the 
LBL appearance measurement of ^13 mentioned in the previous section. We take, based on 
the references, a^ = 0.02/\/2.7 so that S (sin^ 2^13) = 0.02 at 90% CL in the analysis with 
1 degree of freedom. The effective Axl has a definite merit because it does not depend 
on details of the particular experiments, and the value of a^ is a representative of various 
different analyses done before. 
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